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J|h \ Abstract 

The formal extension of the T-duality rules for open strings from Abelian to non-Abelian 
gauge field background leads in a well known manner to the notion of matrix valued 
D-brane position. The application of this concept to the non-Abelian gauge field RG 
/5-function of the corresponding c-model yields a mass term in the gauge field dynamics 
on the matrix D-brane. The direct calculation in a corresponding D-brane model does not 
yield such a mass term, if the Dirichlet boundary condition is implemented as a constraint 
on the integrand in the defining functional integral. However, the mass term arises in the 
direct calculation for a D-brane model with dynamically realized boundary condition. 
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1 Introduction 



T-duality for open strings interchanges Neumann and Dirichlet boundary conditions 
If in a more general setting the string couples with its free ends to an Abelian gauge field, 
the boundary condition stating the balance of the normal derivative of the world sheet 
embedding and the Lorentz force 

d n X» + 2ro'F"„9 t r = 

turns under the T-duality for the i^ 1 coordinate d n X l — ► d t X l (n and t denote normal 
and tangential, respectively) into the condition 

X 1 = 2na'A l {X M ) , 

if the gauge field depends only on a subset X M of the coordinates X^ = (X l ,X M ). 
Thus in the dual description the string endpoints have to move on the manifold, the 
D-brane, defined by the gauge fields in the i-directions. 

The formal extension of this recipe to non-Abelian gauge fields [^] leads to the notion 
of matrix valued D-brane positions which plays a crucial role in the heuristic motivation of 
M(atrix) theory j3J and the emergence of non- commutative geometric structures in string 
theory 0. 

In this paper we will comment on the issue of a cr-model description of the T-dual 
model in the case of non-Abelian gauge fields including its quantum corrections and the 
issue of realizing the boundary condition either as an external one or a dynamical one. 
By external realization we mean a condition on the fields which are integrated over in the 
functional integral, by dynamical realization we have in mind conditions arising as part 
of the stationarity condition of the action. Both kinds of realization are equivalent on the 
classical level, but may lead to different quantum theories. 

The formal manipulations of the functional integral for the cx-model, describing the 
motion of a closed string in nontrivial target space fields (G, B, $) related to its massless 
excitations, which yield Buscher's duality rules ||, have been extended to a a-model 
on a manifold with boundary which in addition couples to a gauge field in ref.0, see 
also ||. There a dual model with externally realized Dirichlet boundary conditions is 
generated. The RG /3-functions of this model have been calculated in ref . |J already. The 
resulting conformal invariance conditions are equivalent to the stationarity conditions of 
the Born-Infeld action dimensionally reduced to the D-brane ||. The Born-Infeld action 
is, up to the dimensional reduction, form invariant under the T-duality transformation 
M. Therefore, at least up to the implemented level of accuracy, the gauge field part of 
the naive T-duality rule is compatible with renormalisation. Q 



With the help of an one- dimensional auxiliary £-field formalism [11], [L2| resolving 



the path ordering prescription for the Wilson loop we were able in M, 13] to extend the 



treatment to the case of non-Abelian gauge fields. In this way the notion of matrix valued 



2 See also the discussion in the second ref. of and the discussion of T-duality and renormalisation 
for the closed string fields JT9 . 
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brane position got a special realization suitable for practical calculations. Postponing 
the ^-integration until the very end, at intermediate stages the non-Abelian gauge field 
appears sandwiched between £ and ( only 

Hence the only difference to the Abelian situation is an explicit dependence on the pa- 
rameter s, parametrising the boundary of the string world sheet. The gauge field and 



brane position RG /3-functions have been calculated in lowest order of a' in [T^ • 

Our only focus in this paper concerns the mass term for non-Abelian gauge field 
dynamics on matrix D-branes which arises by direct application of the naive T-duality 
rule to the /3-function of the dual partner, i.e. the usual open string with free ends moving 
in target space non-Abelian gauge field background. Such a term does not appear in the 
direct calculation of the D-brane model |TB|, [TJJ. Thus at this state of affairs the non- 
Abelian gauge field part of the naive T-duality transformation seems to be incompatible 
with renormalisation at lowest order, already. Before sketching the plan of the paper we 
will recall the formal duality argument giving rise to the mass term 



The lowest order gauge field /3-function for the string with free ends is 0|l5|, [16 



(3{ A) = - a' tT D,F uX . (1) 

If e.g. dxA^ = we get, using X 1 = f(X M ) = -2 7 ra / i 1 (A >A/ ), i] 11 = -1 
and X M = X M , A x = 0, A M = A M , 



^ = -a' rj MN D M F NL + -— - [/, D L f] 

4n z a 



P[ A) = ^V MN D M D N f. (2) 

The second term in turns out as a standard mass term, at least in the special D-brane 
configuration of diagonal and constant /, i.e. /„& = f a S a b, 9m f a = 0, which describes 
planar parallel D-brane copies at position X 1 = f a , a — 1, n. The equation of motion 
obtained from (3^ = then is 

{D M F ML ) ab + (f a - f b f{A L ) ab = . (3) 

As expected the mass is proportional to the separation of the D-brane copies. 

The bulk of the paper is organised as follows. To keep things as clear as possible 
we restrict ourselves to the case of trivial closed string background fields and consider 
T-duality in one direction (A 1 ) only. The calculations crucially depend on the £ auxiliary 
formalism. To become familiar with this calculus we start in section 2 with the derivation 
of the well known variational formulae for the standard Wilson loop |L7| . We do this in 



3 We denote all quantities in this model with a tilde. The untilded quantities refer to the D-brane 
model. 
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a rather explicit way to amplify the step crucially also in the later modified applications: 
The use of the £, £ equation of motion in correlation functions including contact terms. 

To present a simple check of our results in fl3|, [14| free of the subtleties connected 
with the covariant expansion of the imposed boundary constraint, we use in section 3 the 
simplification due to flat target space to integrate the functional 5-function. Things be- 
come even simpler if one then considers constant / only. But also in this very transparent 
situation no mass term appears. 

In section 4 we switch to dynamically realized Dirichlet conditions. The correspond- 
ing model is formulated quantising the action obtained in our paper by a canonical 
transformation of the free end model. As it will turn out, this model yields a gauge 
field /3-function with mass term. In the final section we summarise and interprete our re- 
sult. In particular it will be stressed that the reason for the discrepancy of both versions 
is strongly related to a partial interchange of the functional integration over the string 
position with that over the auxiliary ^-field. 

2 Derivation of the standard variational formula for 
the Wilson loop in (-language 

The Wilson loop as a functional of a closed path X»(s), < s < 1, A^(0) = A^(l) 
defined by 

W[X] = tiPexp(iJ A^X^X" ds^j (4) 

can be expressed by the following functional integral over an one-dimensional auxiliary 
field C, C with propagator (C«(*)C&(s)) = S ab Q(s - t) 0, [12 



W[X] = J D(D( Ca(0)Ca(l) e iSoia] W[X,(,(] ■ (5) 
So and W are defined by 

So[(,c\ = * f'mmds + »c(o)c(o) (6) 



and 

W[X,(,(] = exp^| o 1 CA M (X( S ))C( S )X^ S ) . (7) 

Perhaps it is useful to stress that due to the sum over the index a, realizing the trace, 
there is no dependence on the choice of the point on the closed path from which we count 
s. Varying with respect to X and performing one partial integration we get 

iW (C(<9 M A - d x AJ(X\t) - (A„( - CAX) • 



5Xf(t) 



Now the use of the equation of motion for £ and £ just produces the commutator term, 
necessary to complete the full non-Abelian field strength. To verify this we note that for 
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the full (^-action 

S[(,(] = S [(X]-i\ogW = S [(X]+ f'cA^ds (8) 

Jo 

one has (the corresponding eq. for £ is obvious) 



SS_ 

■5C(s) 

The argument is then completed by the application of 



C = iA^X" - z—^ . (9) 



D(DC C.(0)Ci(l) e i5 l«] h(((t),((t)) T^v = , (10) 

which is valid for an arbitrary polynomial h. Altogether we arrive at [] 

SW 



5X»{t) 



iWCF,xCX\ (11) 



= itrP(e<J>**F M X*(t)) . (12) 

The second functional derivative becomes 

5 2 W - ■ - 

= iW bit - t') Qd u F^C X\t) + iW S(t - t') (F^( (13) 



8X»{t)5X»(t>) "" w 1 "V" "j^w 

- W (F^C X\t) (c(d u A K - d K A„)( X K (t') - (A»( - (A^at: 

In the bilocal (t, i')-term we again want to use the C-equations of motion to get the Abelian 
part of the field strength at if completed by the commutator term to the full non- Abelian 
field strength. But now there is a second factor at parameter value t which is responsible 
for the appearance of a contact term oc 8(t — t') which just covariantizes the derivative 
of F^x in the first term of the r.h.s. of fll3"|). The underlying general formula is (k and h 
polynomials in £, Q 

D(D( C.(0)Ca(l) e iSK '< ] (14) 



/ _ AC _ - r)k \ 







Applying this to eq. (p~3|) we get 

5 2 W 



iW 5(t - f) (DvF^C X (t) + iW 5(t - t') (F^C (15) 



5X»{t)5X v {t') 

- W (F^C X\t) (F VK ( X K {t 



Here and in the following we understand all equations for W and its derivatives up to terms vanishing 
under the £, (^-integration. 
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The transition to sx^)sx"(t') * s ' s i m il ar to the pair (0),(0), obvious and will not be 
written down. Finally we note that (|ITD and (JT3) can be summarised in 

W[X + Y,(,Q = W[X,(X\ (16) 

• exp (i jT 1 ((F, X X X Y» + ^D,F uX X x Y»Y» + ^F^Y^ds + 0(y 3 )) . 

3 /3-functions in the case of external Dirichlet bound- 
ary condition 

We consider the simplest case of only one coordinate of type % and call this X . Then 
the D-brane is defined by X 1 = f(X) with dj = 0. A M (X), diA M (X) = describes 
the gauge field on the brane. Both A and / are matrix-valued. The string world sheet 
M is parametrized by the two-dimensional coordinate z, the boundary dM in z- space is 
described by z(s). The partition function for the model we will discuss in this section is 
given by 

Z[AJ] = J DX»D(D( Ca(0)C«(l) e iSoia] W* 1 " C/POO 

• exp [ d 2 z(dX N 8X N + 3X 1 dX l ) + i [ QA N C, X N ds] . (17) 

\47ra;' Jm JdM J 

It is a special case of the model discussed in ref.|13[. In a first step we can perform the 
^-integration 

i[CK] = J dx 1 exp j M d 2 z OX'dX^j S 9M (X 1 - (f(X)() 

= (det<9 2 H eX p (-?— f X^d^ds 
\47rcr JdM / 



The functional determinant of the Laplacian is taken within the space of functions ap- 
proaching the value zero at the boundary. In the above equation X 1 is the unique solution 
of the boundary value problem 

d 2 X x = in M , X 1 = (f( on dM . (19) 

Up to now there was no need to specify any periodicity condition for £, £. However, since 
the boundary value problem is well posed for C/C(l) = C/C(0) om y> we now choose £ and 
£ anti-periodic |12|]. Then the propagator in the fundamental interval < s < 1 is still 
given by 5 a bQ(s — t). The periodic choice for £, ( would lead to an saw-toothed propagator 
which no longer could organise the path ordering. 

The solution of ([T^) can be represented with some kernel p(z, s) fixed by the geometry 
of M as (For a circle the corresponding integral is the familiar Poisson integral.) 

X\z) = f l dsp(z,s)a S )f(X M (z(s)))((s) . (20) 



3 We switch to Euclidean two-dimensional world volume. 
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Denoting the boundary value of d n p(z,t) by q(s,t) we get (skipping the irrelevant deter- 
minant factor) 

J[C/C] = expf-^/ / (f((s) (fC(t) q(s,t)dsdt) . (21) 
Insertion of this result into (|17]) implies 

Z[AJ] = J DX M D(D( Ca(0)Ca(l) e iSo[c ' fl 

■ exp ^ d 2 z dX 7 ^) • W[X, C, C] , (22) 

where we introduced 

W[XX,(] = ex V (i [ (A N (X N ds - / / (fC(s)(K{t)q(s,t)dsdt) . (23) 

V JdAf 47TO! JdM J<9Af / 

It is very crucial that the modified Wilson loop £, £] is a functional of X (z(s)) 

only. 

Up to now f{X) was an arbitrary matrix- valued function of the coordinates X . To 
keep the variational formulae for W^fX, £, Q as simple as possible, we restrict ourselves in 
the rest of this section to constant f. Then the bilocal term in IV[X, £, Q has influence 
on the (, C equations of motion only. The total action for £, ( instead of S from section 2 
is 

S[C,(] = S [CX}-z\ogW (24) 

= s [c, a + f cA N cx N ds - / I ads) am wsdt , 

Jo 4:7ia' JdM JdM 

leading to 

C(t) = iA N (X N (t)-—K(t) / (f((s)q(t, S )ds - z— . (25) 

lira' Jo <K(s) 

Note that the integral in the last line just represents d n X l . Therefore the equation of 
motion obtained by setting ^|^y to zero coincides with the formal dualization of (|9|). 
Repeating now with the modified quantities the steps presented in section 2 we get 

W[X + Y,(,(] = W[*,C,C] exp(i£l[X,Y,(,(}) , (26) 

with 

o[x,y,c,C] = o(y 3 ) 

* { ((F ML ((s)X L - 2^?C[/, AmK(s) jf 1 (K(t) Y M 

+ ~(D M F NK C(s) X k Y m Y n + ~CF MN ((s) Y M Y N 

-CD M [f, A N ]((s) Y M Y N jf 1 (f((t)q(s, t)dt\ds 



Ana 



47TQ! JO 

1 r ' rl C[/,^]C( s )C[/,^M]C(t)^(^)^ M (t)^,t)M. (27) 



o jo 



From this expression we can read off directly the vertices for a perturbative evaluation 
of Z. Obviously, due to the triviality of our target space fields (G, B, there are no 
vertices in the bulk of the string world sheet, but on the boundary only. The propagator 
(Y M (s)Y N (t)) restricted to the boundary is equal to — 2a'rj MN log \s — t\. Since this is 
integrable, in lowest order the only divergences are due to the tadpole diagrams arising by 
connecting the y-legs of either the third or the fifth vertex in ( p7D . The divergence arising 
from the third term contains X , hence it contributes to the gauge field /^-function. The 
fifth term contains dnX 1 = j (f(q(s,t)dt. This is an independent divergence |9], [L3| and 
constitutes the /3-function for the brane position /. Altogether we find 

= -a'D M F MN 

P if) = ^D M [f,A M ] = ±D M D M f. (28) 

4 /^-functions in the case of dynamically realized 
Dirichlet boundary condition 

In this section we consider a model defined by a functional integral without any constraint 
on the functional integrand and take as the action that which arises by formulating the 
T-duality as a canonical transformation 0. Introducing 

W[X,(,Q = expij dM ((A N (X N + -^(Kd Il X^ds, (29) 

the partition function is given by 

Z[AJ] = J DX P D(D( Co(0)Co(l) e iSoia] 

■ exp (-^ f d 2 z dX^dX^ ~- - / X'd^X.ds) W[X, C, C] • (30) 
\47rcr Jm 2ira' JdM J 

The Dirichlet condition arises as part of the stationarity condition of the action since 
X 1 — (f( is the factor multiplying d-^dXi on DM. In contrast to the Wilson loops of the 
previous sections W is, due to the presence of d n Xi, a functional of X^(z) not only on 
the boundary, but on the string world sheet M itself. Again, after a partial integration, 
we get as a starting point for the variational formulae 

rtU i W C 5 {2 \z(s) -z) (31) 



5X N {z) Jo 



C(d N A K - d K A N )(X K + ^i&NfCdnX! - (A N ( - (A N ( j d, . 



The total action for £, ( is 



%,C] = S [(,(]-ilogW (32) 

= S [(, C] + C (A N (X N ds + / (Kd^ds . 
Jo 2na' JdM 
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With 



i <W 
((t) = iA N (X N (t) + —fffldnX! - (33) 



we repeat the procedure of section 2 and get 



SW 

5X N (z) 



iW £ 6W{z{s) - z) ((F NK C X K + ^(DnK d n X^ ds , (34) 



• f 1 5®(z(t) - z') ((F ML (X L + J- 1 CD M fCd n X 1 ) dt 
Jo \ 2na J 

+iW J 5 (2 \z(s) - z')5 (2 \z(s) - z) ((D M F NK (X K + ^— t (D M D N f(d n X^) ds 

+iW [ (F NM ({s) 6®(z(s) - z) 4-5® {z(s) - z')ds . (35) 
J ds 

The variations with respect to X 1 are different qualitatively. There is no way to get rid 
of d n acting on the variation of X 1 by some partial integration on the boundary. Thus 
no C or ( is generated, and we have no opportunity to apply the equations of motion for 
C, C i n a useful manner. | The first and second derivatives are 



SW lW - jlfidjl- 2 \z(s)-z)ds, (36) 



5Xx{z) 2vra' 
5 2 W W 



There is still a mixed derivative sx 1 (z)sx N (z') ' we save wr iting down its form and 
proceed to the combination of all first and second derivatives in 

W[X + Y, C, C] = W[X, C, C] exp(zfi[X, Y, C, CD , (37) 

with 

n[x,Y,c,C) = J o 1 {(CF N LC{s)x L + ^ ; CD N fCdnX 1 )Y 1 

+ ?L>~ uc dnYi + ^ Fnm<: fMyN 

+^((D M F NK ( X K + -L-(D M D N f( d^X^Y 
+^CD M f( Y M d n Y 1 + 0{Y 3 )}ds . (38) 



N 



5 Fortunately here these manipulations are not needed to get manifest gauge invariant structures. 



s 



To calculate the partition function ( p0|) we make a shift X — > X + Y, use -D(X + y) = DY 
and choose X in such a way that the linear terms in Y vanish in the bulk of M, i.e. 
d 2 X fl = 0. However, we do not specify any boundary condition on X^. [] Then we have 
to continue with 

Z = |dODC Ca(0Ra(l) e iSo ^ + ^f( xN9 " XN - xld » x ^ ds W[X,(,(] (39) 
■ J DY V exp (J dY^dY^z -2 J Y 1 d a Y l ds^j\ 

•expj^ J (Y]\fd n X N — X 1 d n Yi)ds + tQ[X,Y,C,C]] ■ 
For the perturbative evaluation of the F-integral it is convenient to use 

= (detA)5exp|^ J J a (z)A al3 (z,z , )J^(z')d 2 zd 2 z^ . (40) 

The propagator A of the above equation is defined by [] 

A^ = ?7 Mi , A u (no sum) , 

d 2 A^ v = -2<Ka l %J { - 2) (z- z') , inM, A x = ,d n A N = , on dM . (41) 

Eqs. (|39D-(|4lD imply a set of Feynman rules with Neumann propagator for Y N , Dirichlet 
propagator for Y 1 , no vertex in the bulk and boundary vertices defined by Cl as well as 
the two additional 1-leg vertices in the last exponential of eq. . 

The two vertices in the third line of eq. fl38|) yield divergent tadpole contributions which 
would deliver just the same /3-functions as in section 3, if they would be the only divergent 
diagrams. But in contrast to section 3, where 1-leg vertices appeared with underivated 
Y only, now we have the 1-leg vertex 2^y(C/C — X l )d n Yi. The diagram connecting two 
such 1-leg vertices is 

A = \ (2^p / ^ ( ^ /C(s) " Xl )GK® ~ Xl) d - d 'n A ^ *')U> ■ ( 42 ) 
For the upper half plane 

Ax{z,z') = - a' (\og\z - z'\ - log|z-?|) (43) 

leads to 

-i f (C/CW - *')(</««) - x 1 



a = ^i dsdt (s _y — ■ (44) 



7 We could choose X in such a way that all linear terms in Y would vanish, but then the factorizing 
classical factor would contain £, £ and mix with W[X,(,C] an( i the Y functional integral in doing the 
(^-integration. 

8 To avoid the discussion of the zero mode for the Neumann propagator we take M to be the upper 
half plane. 
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A is linearly divergent. We neglect linear divergences f\ and look only for logarithmic ones 
which could be produced by expanding the nominator in A. If this nominator would be 
continuously differentiable, due to the antisymmetry of (s — £) -1 , there would be no such 
divergence. However, we have to evaluate the diagram with quantised fields Y and £, (. 
As a consequence the left and right limits of the derivative of the nominator could differ. 
Anticipating this possibility we expand the nominator for s > t and s < t separately 



C/COOC/CW = 0(1) + 0»-f){e(a-f) ^ (C/C(* + o)) C/CW 

+ e(t - a) I (C/C(* - 0)) CK(t)} + 0((a - t) 2 ) . (45) 
Using the equations of motion ( |33f) this turns into 

C/COOC/CC*) = 0(i) + («-t)^{e(a-t)C^/C(*+o)C/C(*) 

+ e(t - fl ) c^iv/C(t - o) c/C(*)} + o(( s - 1) 2 ) . (46) 

Now from the special structure of the ^-propagator it is obvious that for the purpose 
of the perturbative evaluation of the £> C functional integral we can make the following 
identifications 

C«(* + o) C&(* + o) Cc(*) Cd(t) = he L{t + o) Cd(t) , 

Ut - 0) C*(* - 0) Cc(*) Cd{t) = S ad Cc{t) Ut - 0) . (47) 



Using this in (pE6| ) we get for the regularised version (a regul. parameter) of fl4~4| ) Q 



?<t (s — t) 2 + a 2 

the dots standing for either linear divergent or finite contributions. This implies 

A = ^log« / X N ([f,D N f]((t)dt + ... . (48) 

This divergence is responsible for the presence of the mass term in the gauge field /3- 
function. Altogether we find 

^ = -a' D M F MN + ^[f,D N f] . (49) 

coincides with that of section 3. 



9 They are irrelevant in dimensional regularisation. 

10 The terms in the nominator of ( f4~i| ) containing f in zero or first order are continuously differentiable. 
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5 Conclusions 



At least up to the considered order of perturbation theory, our results in section 3 and 
4 clearly favour the D-brane model with dynamically implemented Dirichlet boundary 
conditions as the correct realization of the T-dual partner of the open string with free 
ends. The model with externally realized boundary conditions has been motivated |7], [13 



by formal manipulations of the functional integral. Its failure to produce the mass term 
can be explained just by this formal nature. However, on a more ambitious level we 
would like to localise the forbidden step in those manipulations. A first idea we get by 
looking back to the application of the naive T-duality rules to the gauge field /5-function 
of the open string with free ends. The mass term arises from taking the value 1 for 
the summation index. The underlying ultraviolet divergence is due to a propagator in 
1-direction at coinciding points. Since there is no propagation in this direction in the 
model of section 3, it should be not surprising that no mass term appears. On the more 
technical level the comparison of the details of the calculations is useful. 

In both cases we consider the partition functions Z^ Z as defined by the interacting 
quantum field theory of the string world sheet position X and the auxiliary fields £, ( in 
the sense of a combined functional integral. This point of view is implemented by the use 
of the ( equations of motion including contact terms for the expansion of the modified 
Wilson loops. However, we have given up this understanding at one occasion, namely 
when we switched to a partially stepwise evaluation by performing the X 1 integration in 
section 3. But just the intermediate treatment of £, ( as classical fields is the dangerous 
step as can be seen in the calculation of diagram A in section 4. If we there had treated 
C, C as classical fields we would also have lost the mass term. 

The stepwise, i.e. first X and then £, £ integration has been crucial for giving the 
formal notion of matrix valued brane position a technical meaning via insertion of matrices 
in the argument which specifies an explicit boundary parameter dependent scalar valued 
Dirichlet boundary condition |7|, [13|. Now we see that just this stepwise procedure violates 
the equivalence with the original free end model. Nevertheless the model of |TB[] is well 
defined and interesting for its own. 

We stress that the preference of the dynamically over the external realization of bound- 
ary conditions developed in this paper concerns the case of non-Abelian gauge fields and 
the implementation of path ordering via the ^-formalism only. In the Abelian case we see 
no stumbling block in the equivalent transformation of the free end model to the D-brane 
model with externally realized boundary conditions. We can also not exclude that in 
the non-Abelian case there could exist T-dual versions with external boundary conditions 
beyond the ^-formalism, for some alternative aspects of the construction of a cr-model 
describing matrix D-branes see |I8|. 



Our discussion concerned the test of quantum equivalence of naive T-dual models 
in the general, not necessary conformally invariant situation. Since T-duality as a map 
between equivalent theories is by far better established in the conformal case, it would 
be interesting to understand both the potential differences between externally and dy- 



ii 



Note that in both cases the £ equations of motion are compatible with naive duality. 
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namically realized boundary conditions as well as the a-model description of boundaries 
coupling to non-Abelian gauge fields in the formalism of boundary conformal field theo- 
ries, and refs. therein. 

Acknowledgement: I would like to thank H.-J. Otto, C. Preitschopf, A. Recknagel and 
V. Schomerus for useful discussions. 
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